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Abstract 

Exclusive semileptonic B meson decays to orbitally excited D mesons are 
investigated in the infinitely heavy quark limit in the framework of the rela- 
tivistic quark model based on the quasipotential approach. The B — > D**ev 
Isgur-Wise functions t 3 / 2 ( w ) an d Ti/2( w ) are determined. It is found that 
the relativistic transformation of the meson wave functions (Wigner rotation 
of the light quark spin) contribute already at the leading order of the heavy 
quark expansion. 

The investigation of semileptonic decays of B mesons to excited D meson states is a 
problem which is important both from a theoretical and experimental point of view. In 
particular, these decays can provide an additional source of information for the determination 
of the Cabibbo-Kobayashi-Maskawa matrix element V c b as well as on the relativistic quark 
dynamics inside heavy-light meson. The experimental data on these decays are becoming 
available now 10-0], and the B factories will provide more accurate and comprehensive 
data. The presence of the heavy quark in the initial and final meson state in these decays 
considerably simplifies their theoretical description. A good starting point in this analysis 
is the infinitely heavy quark limit, itlq — > oo [f|]. In this limit the heavy quark symmetry 
arises, which strongly reduces the number of independent weak decay form factors |J. The 
heavy quark mass and spin decouple then and all meson properties are determined by light 
degrees of freedom alone. As a result the heavy quark degeneracy of levels emerges. The 
spin s q of the light quark couples with its orbital momentum I (j — I ± s q ), resulting for 
P-wave mesons in two degenerate j = 3/2 states ( J p = 1 + , 2 + )[] and two degenerate j = 1/2 
states (0 + , 1 + ). The heavy quark symmetry also predicts that the weak decay form factors 
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1 Here J = j ± 1/2 and P are the total angular momentum and parity of the meson. 
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for B — > D**ev decays, where P** is a generic P-wave D meson state, can be expressed in 
terms of two independent Isgur-Wise functions r 3 / 2 and T\/2 @- 

In preceding papers we have calculated the mass spectra of orbitally and radially excited 
states of heavy-light mesons as well as different weak decays of B mesons to ground 
state heavy and light mesons ]7|,|8| in the framework of the relativistic quark model based 
on the quasipotential approach in quantum field theory. Let us now apply this model to the 
investigation of semileptonic B decays to P-wave D mesons in the heavy quark limit. 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation || of the Schrodinger 
type [0 

where the relativistic reduced mass is 

"* = Wr- 91 ' (2) 

Here m qt Q are the masses of light and heavy quarks, and p is their relative momentum. In 
the centre-of-mass system the relative momentum squared on mass shell reads 

h % (M , [M 2 - (m q + m Q f] [M 2 - (m q - m Q ) 2 } 

b{M) ~ AM 2 • {6) 

The kernel V(p, q; M) in Eq. (|I]) is the quasipotential operator of the quark- ant iquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by the 
Lorentz-structure of the confining quark-antiquark interaction in the meson. In constructing 
the quasipotential of quark-antiquark interaction we have assumed that the effective inter- 
action is the sum of the usual one-gluon exchange term and the mixture of vector and scalar 
linear confining potentials. The quasipotential is then defined by fil l 



V(p,q; M)=u q (p)u Q (-p)^a s D flu (k)^Y Q 

+v;L(k)r^r Q;M + Kl f (k) j^(g) MQ (-g), (4) 

where a s is the QCD coupling constant, D^ v is the gluon propagator in the Coulomb gauge 
and k = p — q; 7^ and u(p) are the Dirac matrices and spinors 



u x {p) 



e(p) + m 




(5) 



- t{p)+m 



with e(p) = Vp 2 + m' 2 . The effective long-range vector vertex is given by 



r„(k) = 7 M + 2^<v r > ( 6 ) 
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where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

= (1 - e)Ar, 

V c s oni (r)=eAr + B, (7) 

reproducing 

Konf (r) = Kif (r) + Konf (r) = At + B, (8) 

where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, expanded in v 2 /c 2 , can be found in 
Refs. Q,Q and for heavy-light mesons in |J. All the parameters of our model like quark 
masses, parameters of the linear confining potential, mixing coefficient e and anomalous chro- 



momagnetic quark moment k were fixed from the analysis of heavy quarkonia masses [11 



and radiative decays [O]. The quark masses m& = 4.88 GeV, m c = 1.55 GeV, m s = 0.50 
GeV, m u ^ = 0.33 GeV and the parameters of the linear potential A = 0.18 GeV 2 and 
B = —0.30 GeV have usual quark model values. The value of the vector-scalar mixing coef- 



ficient £ = — 1 has been determined from the consideration of the heavy quark expansion [14 



and meson radiative decays [ 13| . Finally, the universal Pauli interaction constant k = — 1 



has been fixed from the analysis of the fine splitting of heavy quarkonia 3 Pj- states [|TT] . 
Note that the long-range magnetic contribution to the potential in our model is proportional 
to (1 + k) and thus vanishes for the chosen value of k — — 1. 

In order to calculate the exclusive semileptonic decay rate of the B meson it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
The matrix element of the weak current J w = C7 M (1 — 7 5 )6 between B meson and orbitally 
excited D** meson in the quasipotential method has the form [JTSJ 

(b~\J?(o)\b) = J $^*Mp)r M (p,q)**(q), (9) 

where r M (p, q) is the two-particle vertex function and ^b,d** are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to T come from Figs. 1 and 2.0 In the heavy quark limit rrif, c — > oo 
only contributes, while contributes at 1 /iiiq order. As we limit our analysis here to 
the leading order of the heavy quark expansion, only the vertex function is necessary. 
It looks like 

r«(p, q) = u c ( Pc )l,(l - 7 5 H(g & )(27r) 3 5(p g - q q ), (10) 

where UI5 



2 The contribution I^ 2 ) is the consequence of the projection onto the positive-energy states. Note 
that the form of the relativistic corrections resulting from the vertex function I^ 2 ) is explicitly 
dependent on the Lorentz-structure of the gg-inter action. 
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Peg = e c , q (j>)j^±J2n®(p D **)p i 

D** i = i 

q b , q = %q(p)jj- ± E nW (Pfi)?\ 



i=l 

and 



M' J M(E + M)j ' 
The wave function of a P-wave P** meson at rest is given by 

*^(p) = *^,(p) = y/ w ^ w (p) ) (ii) 

where J and M are the total meson angular momentum and its projection, while j is the 
light quark angular momentum. ^du){p) is the radial part of the wave function, which has 
been determined by the numerical solution of eq. ([]]) in [0. The spin-angular momentum 
part has the following form 

y J 3 M = J2UM-a Q , la Q \JM)(lM-a Q -a qi \ a q \j M - a Q )Y^-^ X Q^Q)xM- 

(12) 

Here (jimi, j2"m2\J M) are Clebsch-Gordan coefficients, Y™ are spherical harmonics, and 
xip) (where a = ±1/2) are spin wave functions, 

x(l/2)= Q, Xi-l/2) = (^j. 

It is important to note that the wave functions entering the weak current matrix element 
(§) are not at rest in general. E.g., in the B meson rest frame the P** meson is moving with 
the recoil momentum A. The wave function of the moving P** meson ^> £>** a is connected 
with the P** wave function at rest ^d**o = ^D(j) by the transformation [15| 

a(p) = ^Wa^Wa^-oCp), (13) 

where R w is the Wigner rotation and the rotation matrix P 1//2 (P) in spinor representation 
is given by 

(J °^Dl/ q 2 (RY A ) = S-\ Pc , q )S(A)S(p), (14) 

where 



<?(p) = >(P)+™ ( 1 + °P ) 
V 2m y e(p) + m J 

is the usual Lorentz transformation matrix of the four-spinor. For electro-weak B meson 
decays to >S-wave final mesons such transformation contributes at first order of the 1/mq 
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expansion, while for the decays to excited final mesons it gives a contribution already to the 
leading term, due to the orthogonality of the initial and final meson wave functions. 

In the infinitely heavy quark limit (m& jC — > oo) all form factors of the semileptonic 
B — > D**ev decays are related to two independent Isgur-Wise form factors T\/ 2 an d t 3 / 2 by 



(D (l/2)(tO|&y M (l - l 5 )b\B(v)) = 2 v /M B M Do(1/2) r 1/2 (w)(v'^ - v,), 
(D 1 {l/2){v')\c lfi {l - ! 5 )b\B{v)) = 2^M B M Dl{1/2) T l/2 {w){ie mPl t* a v p v^ 

+ (1 - wK + (e* • 

(D a (3/2)(t0i&y M (i - i 5 )b\B(v)) = 



+(l-^) e *-( e *. w )[3^-( w -2W]}, 



(D 2 (3/2)(0|c 7/ ,(l - l 5 )b\B(v)) = ^3M B M D2{3/2) r 3/2 (w){ - i^e^^A'^ 



(15) 



where f (v') is the four-velocity of the initial (final) meson, w = v ■ v', and e M , are 
polarization vector and tensor of D\ and D 2 mesons, respectively. 

To calculate the corresponding matrix elements, we substitute the vertex function 
( [TOD in the matrix element of the weak current between meson states ([]) and take into 
account the wave function properties (|TT|)-(|i"3|). Then, in the limit m^^ — > oo we find that 
the heavy quark symmetry relations (|15|) are exactly satisfied in our model. The resulting 
expressions for the Isgur-Wise functions r 3 / 2 and Ty 2 are 



T3/2IWJ 



V2 



1 



x 



3 (w + iyi 2 J (2tt) 3 
d p 



d 3 p - 

^D(3/2)(P 



2e„ 



-2e. 



Ti/ 2 (w) 



1 



9 dp 
1 



+ 



^b(p), 



Sx/^^ + l) 1 / 2 J (2tt)3 
2p 



d 3 p - 

VJJ(1/2)(P 



Mfl( 3 /2)(«) + 1) 



2e, 

M D(l/2) (w + 1)' 



(16) 



x 



-2 — 

q dp 



e q + m g 



^b(p), 



(17) 



where the arrow over 8/ dp indicates that the derivative acts on the wave function of the 
D** meson. The last terms in the square brackets of these expressions result from the wave 
function transformation (|T3|) associated with the relativistic rotation of the light quark spin 
(Wigner rotation) in passing to the moving reference frame. These terms are numerically 
important and lead to the suppression of the T1/2 form factor compared to r 3 / 2 . Note that if 
we had applied a simplified nonrelativistic quark model |5,16] these important contributions 
would be missing. Neglecting further the small difference between the wave functions ipD(i/2) 



and ipD(3/2), the following relation between r 3 / 2 and T\/ 2 would have been obtained |T7 



n/ 2 (w) 



w + 1 



-T3/2{W) 



(18) 
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However, we see that this relation is violated if relativistic transformation properties of wave 
function are taken into account. At the point w = 1, where the initial B meson and final 
D** are at rest, we find instead the relation 



r 3/2 (l) - r 1/2 (l) - ~ / tI^Wp)-^— ^(P), (19) 



2 J (2tt) 3Y ^ ^ J e q + m q 

obtained by assuming ipD(3/2) — i>D(i/2) — ^d**- This relation ( p79|) coincides with the one 
found in Ref. [TJ| where the Wigner rotation was also taken into account. 



In Table | we present our numerical results for 7^(1) and its slope p 2 = — ~^7j 



m 

UI = 1 



comparison with other model predictions |T7|-p3|. Moreover, we plot tj(w) for B — ► D**ev 



and B s — ► D**ev decays as function of w in Figs. 3, 4. The corresponding decay rates and 
branching ratios are given in Tables || and |HI[ We see that most of the above approaches 
predict close values for the function r 3 / 2 (l) and its slope p^ 2 , while the results for Ti/ 2 (1) 
significantly differ from each other. This difference is a consequence of a different treatment 
of the relativistic quark dynamics. Nonrelativistic approaches predict t 3 / 2 (1) — 71/2(1) (see 
(fjjD), while the relativistic treatment leads to 73/2(1) > 7"i/ 2 (l) (see (|T9|)) . Our results for 
the branching ratios of B — > Di j2 (3/2)ez/ decays are consistent with available experimental 
data which at present require to use some assumptions about the branching fractions 
of the Dj mesons. 

Finally, let us test the fulfilment of the Bjorken sum rule |M] in our model. This sum 
rule states 



P 2 = \ + H \r$(l)\ 2 + 2£ |r^(l)| 2 + ■ ■ • , (20) 



ni 



where p 2 is the slope of the B — > D^eu Isgur-Wise function, Ty£ an d T 3/2 are ^ ne f° rm 
factors describing the orbitally excited states discussed here and their radial excitations, and 
ellipses denote contributions from non-resonant channels. We see that the contribution of 
the lowest lying P-wave states (m = 0) implies the bound 

P 2 > ~ + ki/ 2 (l)| 2 + 2|r 3/2 (l)| 2 = 0.80, (21) 



which is in agreement with the slope p 2 = 1.02 in our model [14 ]. 



In this paper we have applied the relativistic quark model to the consideration of semilep- 
tonic decays of B mesons to orbitally excited charmed mesons in the leading order of the 
heavy quark expansion. In particular, it has been found that the Lorentz properties and 
transformations of meson wave functions play an important role in the theoretical descrip- 
tion of these decays. Thus, the Wigner rotation of the light quark spin gives a significant 
contribution already at the leading order of the heavy quark expansion. In conclusion let us 
mention that the corrections in inverse powers of the heavy quark mass l/m c ^ to the decay 
rates might turn out to be non- negligible, especially for spin zero and spin one D** mesons 



17| . The investigation of such corrections in the framework of our model is an important 
task that will be considered elsewhere. 
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TABLES 



TABLE I. The comparison of our model results for the values of the functions Tj at zero recoil 
of final D** meson and their slopes p\ with other predictions. 



our 







73/2(1) 
P3/2 

n/ 2 (i) 

Pl/2 



0.49 
1.53 

0.28 
1.04 



0.41 
1.5 
0.41 
1.0 



0.56 
2.3 

0.09 
1.1 



0.66 
1.9 
0.41 
1.4 



0.35 ± 0.0* 
2.5 ±1.0 



0.54 
1.5 
0.22 
0.83 



0.52 
1.45 
0.06 
0.73 





TABLE II. Decay rates L (in 


units 


| Kb/0.041 2 10~ 15 GeV) and branchhn 


I ratios (%) for 


B 


— > D**ev decays. 












Decay 


L 


Br 


Br (CLEO) § 


Br (ALEPH) § 


B 


-> Di(S/2)eu 


1.4 


0.33 


0.56 ±0.13 ±0.08 ±0.04 


0.74 ±0.16 


B 


-> D 2 (3/2)ev 


2.1 


0.52 


< 0.8 


< 0.2 


B 


-►Z>i(l/2)ei/ 


0.30 


0.074 






B 


A)(l/2)ei/ 


0.25 


0.062 









TABLE III. Decay rates L (in units 


Vcfe/0.04 2 10 15 GeV) and branching ratios 


;%) for 


B s 


— > D**eu decays. 








Decay 


r 


Br 


B s 


-► Aa(3/2)e*/ 


1.6 


0.39 


B s 


D s2 (3/2)e*/ 


2.5 


0.59 


B s 


-► £> s l(l/2)ei/ 


0.54 


0.13 


B s 


-> £> s o(l/2)e!/ 


0.45 


0.11 
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FIGURES 




FIG. 1. Lowest order vertex function contributing to the current matrix element (9). 
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FIG. 2. Vertex function r^ 2 ^ with the account of the quark interaction. Dashed lines correspond 
to the effective potential (||). Bold lines denote the negative-energy part of the quark propagator. 
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FIG. 3. Isgur-Wise functions t s / 2 (w) (upper curve) and Ti/ 2 (w) (lower curve) for B — ► D**ev 
decay. 




